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Abstract 

Due to a first order phase transition, a compact star may have a discontinuous distribution of 
baryon as well as electric charge densities, as e.g. at the surface of a strange quark star. The 
induced separation of positive and negative charges may lead to generation of supercritical electric 
fields in the vicinity of such a discontinuity. We study this effect within a relativistic Thomas- 
Fermi approximation and demonstrate that the strength of the electric field depends strongly on 
the degree of sharpness of the surface. The influence of strong electric fields on the stability of 
compact stars is discussed. It is demonstrated that stable configurations appear only when the 
counter-pressure of degenerate fermions is taken into consideration. 

PACS numbers: 97.60.-s, 12.38.Mh, 04.70-s 
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Physics of compact stars is under intensive investigation for many decades. Modern 
studies are focused on the equation of state (EOS) of dense baryonic matter including 
possible phase transitions to the quark matter. One interesting feature of a strong first 
order phase transition is that it can generate a sharp discontinuity in the baryon density as 
a function of radial coordinate, see e.g. ref. [1]. In the extreme situation when the EOS 
has a zero pressure point at a finite baryon density pc, as is the case e.g. for the MIT bag 
and NJL models [2, 3], the discontinuity in the baryon density occurs at the surface of the 
star, i.e. p{R) — Pc- Such a situation is expected in compact stars made of Strange Quark 
Matter (SQM), as first predicted in refs. [4, 5]. If the matter is composed of several species 
with opposite electric charges and different masses, the presence of a sharp discontinuity 
should lead to a charge separation and generation of an electric field. This effect was first 
discussed in ref. [6] in context of the necked SQM stars, where the electrons from the bulk 
SQM matter can penetrate in vacuum through the sharp star's surface. As demonstrated 
in ref. [7], supercritical electric fields can be generated at the boundary of a nuclear core 
surrounded by an electronic cloud. In ref. [8] the structure of compact stars with a net 
charge at the surface was calculated within the General Relativity. The boundary effect may 
be important also for understanding the structure of a mixed phase, where domains of two 
phases have opposite electric charges [9, 10]. The electrostatic interactions are important for 
the description of neutron-star crusts where atomic nuclei are embedded in a dense electron 



Our goal in the present paper is to study the conditions leading to the generation of 
supercritical fields at a star boundary. We solve a simplified problem replacing the spherical 
star boundary by the planar one. We expect this approximation to work very well for a 
compact star with radius of about 10 km. We assume further that the net positive charge, 
associated with protons or quarks, has a smooth-step distribution of the Woods-Saxon type 



where is the positive charge density in the bulk [z — j- — oo) and a is a diffuseness 
parameter. The star boundary is located at z — zq and Pp{zq) — Ppo/2. At a ^ this 
distribution approaches a rectangular step considered in ref. [7]. 

For a macroscopic object hke a star the condition of global charge neutrahty must hold 
to a very high precision, see discussion in ref. [13]. Therefore, the positive charge of Eq. 



gas [11, 12]. 
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(1) must be fully neutralized by the negative charge of electrons. However, since electrons 
are light and interact only via the electromagnetic force, they will penetrate through the 
boundary and generate a local charge disbalance around the star surface. The induced 
electrostatic potential (j){z) is determined from the Poisson equation (see e.g. [14]) 

(P(f) 

— = -e [pp{z) - peiz)] = -epch(^) , (2) 

where e = -\/47ra=0.3028 is the proton charge and a is the fine-structure constant[26]. 
For a given proton distribution, Eq. (1), the electron charge distribution Pe{z) should be 
determined self-consistently. For this purpose we use the Thomas-Fermi approximation 

[15, 16], which should work well for an extended object like a heavy nucleus or star. The 
relativistic version of this method was considered e.g. in rcfs. [17, 18]. In a semi-classical 
approximation the electron energy at point z can be written as 

e{k,z) = ^yk^ + ml-V{z) , (3) 

where k is its 3-momentum and —V{z) = —e(f){z) is the potential energy. At zero tem- 
perature all electronic states with k < kpi^z) are occupied, where kp^z) is the local Fermi 
momentum. It is determined from the condition 

eikp^z), z) = \J kp{z) + ml — V{z) = fi — const , (4) 

that gives 

kF{z) = ^[^ + V{z)f-ml . (5) 
The local electron density is found by integration over k: 



, , ^ f d'k [{p + V{z)y-mll 
P^^'^^^ J (2^^ 3^^ 





(6) 



By inserting this expression into Eq. (2) we obtain a non-linear differential equation for 
To find its solutions wc need to specify boundary conditions. 
According to our assumption, the boundary is in direct contact with the vacuum and 
therefore the electron density as well as the electric potential must vanish at 2; ^ oo. This 
condition can be fulfilled if // < mg. For p < rrie the electrons are bound to the surface, i. e. 
a finite energy is needed to extract an electron from the star. This energy is well known in 
ordinary metals as the exit work. In this case the electron density should vanish at a certain 



distance from the surface, z*, which is found from the condition nie — V{z^) — /i, see Eq. 
(4) . The value of /i is determined from the condition of global charge neutrality: 

oo 

e J [Pp{z)-pe{z)]dz = a+ + a^^O, (7) 

— oo 

where cr+ and (T_ are the areal densities of positive and negative charges accumulated at 
z < Zq and z > Zq respectively. Since calculated /i values are very close to mg, we use /i — rrie 
for our estimates. 

Another boundary condition follows from the requirement that far inside the star the 
electron charge density must completely compensate the positive charge density, i.e. Pch{z — > 

1 /3 

— oo) = 0, that means that kpei^ -> —oo) = kpo = (Svr^ppo) ^-nd 

V{z ^ -oo) = ^^1^0 + ml -11^ ^(37r%o)'/' + m2 - . (8) 

These conditions are sufficient to determine the electron density distribution Pe{z). Unfor- 
tunately, no analytical solution can be found without further approximations, and therefore 
numerical methods must be applied. 

We have performed numerical calculations for different values of the diffusencss parameter 
a ranging from 1 fm, typical for nuclear surfaces, up to 1 nm, typical for solid-state surfaces. 
As an example, in Fig. 1 we present results for a—h fm which is rather close to a "sharp 
surface" limit. The upper, middle and lower panels show, respectively, the electron density, 
electric potential and electric field calculated for the given proton density (dashed hue in the 
upper panel). One can clearly see the deviation of the electron density from the proton one in 
the layer of ±50 fm around the surface. This gives rise to a sharp fall of the potential, and as a 
consequence, to a very strong electric field generated in this layer. The maximum of the field, 
about 1.5 MV/fm, is reached Bi z — Zq. This value is almost 600 times higher than the critical 
electric field for spontaneous electron-positron pair production [19], Ec — 2w?c^ /eh — 0.0026 
MV/fm=2.6 • 10^^ V/cm. However, in the considered system the pair production does not 
happen because all electron states with ee{k,z) < rUe are occupied and therefore Pauli- 
blocked. 

This situation is clarified in Fig. 2, where the electron levels in a strong electric potential, 
—V{z), are shown schematically. As explained above, the negative potential inside the 
medium {z < zq) is generated by the charge separation at the surface. In the vacuum such 
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FIG. 1: The upper, middle and lower panels show, respectively, the electron density (full line), 
electric potential and electric field, calculated for the given proton density profile (dashed line in 
the upper panel) as given by Eq.(l) with ppQ = 0.01 fm~^ and a=5 fm. 

a potential would certainly lead to a spontaneous electron-positron pair production. In 
this process the negative-energy electrons from the Dirac sea (below the mass gap) would 
penetrate through the barrier and occupy a state in the upper continuum (above the mass 
gap) leaving a hole (positron) behind. But all electron states up to e = are occupied 
(Fermi sea) to balance the bulk positive charge of protons. The situation is changed at 
finite temperatures when the particle-hole excitations are produced in the Fermi sea. Then, 
the electrons from the Dirac sea can occupy the vacancies in the Fermi sea. The process of 
thermal emission of electron-positron pairs from the sharp surface was considered earlier in 
ref. [20]. Here we want to point out other interesting processes, namely, the spontaneous 
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FIG. 2: Schematic view of the electron energy levels in a strong electrostatic potential —V{z). The 
electron states above (Fermi sea) and below (Dirac sea) the mass gap (white area) are shown by 
dark and grey shadowing. Transitions from the Dirac to Fermi sea (long arrow) are not allowed 
because all states in the Fermi sea are occupied. 

production of negative muons and pions, which may replace the electrons. They become 
possible in a strong enough field when Vq > mfj_, rriT^. This opens the possibility of pion 
condensation induced by strong electric fields [17]. 

With increasing a the electron density comes closer and closer to the proton density 
and the charge separation at the surface diminishes. The electric potential difference at 
z — )■ ±00 is fixed by the boundary conditions, but its gradient, i. e. the electric field, across 
the transition region becomes smaller and smaller. We have checked that the electric field 
falls below the critical value at a ~ 10 pm. In Table 1 we give the maximum value of the 
electric field, Emax, calculated for several values of the diffuseness parameter a. Emax can be 
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TABLE I: Maximum electric field and areal charge density calculated for several values of the 
diffuseness parameter a. 



a 


E 

'-'max 


[MV/f„,] 




<r+ [e/f 




1 fm 


1.82 


•10-0 


6.89 


10+2 


1.01 


• 10" 


-1 


10 fm 


1.16 


•10-0 


4.38 


10+2 


6.41 


• 10" 


-2 


100 fm 


2.05 


•10-1 


7.74 


10+1 


1.13 


• 10" 


-2 


1 pm 


2.08 


•10-2 


7.86 


10 


1.18 


• 10- 


-3 


10 pm 


2.08 


• 10-3 


7.86 


10-1 


1.18 


• 10- 


-4 


100 pm 


2.08 


• 10-^ 


7.86 


10-2 


1.18 


• 10" 


-5 


1 nm 


2.08 


• 10-5 


7.87 


10-3 


1.25 


• 10- 
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calculated by using the formula which follows from the Gauss law: 



Emax = E(2; — Zq) — e 



zo 

J Pch{z)dz 



(9) 



where a+ is the areal density of the net positive charge at ^ < 

For our discussion below we need to estimate the electric field and other related quantities 
on a qualitative level. This can be done by using the ultrarelativistic limit, V ^ me, which 
is justified everywhere except of the asymptotic region at z — )■ oo. Taking for simplicity the 
hmit a — > 0, we get the exact solution of the Poisson equation (2) [27] 



V{z) 



l-ciexp(^^ 



C2 + 



V6A 



Z < Zo 
Z > Zo 



(10) 



where 



A- 



71 



ek 



FO 



(11) 



The constants ci = 0.2374 and C2 = 1/(1 — ci) = 1.3113 are obtained from the continuity 
conditions for V{z) and dV/dz at z — zo- It is interesting that the decay length A is 
determined by the electron Fermi momentum in the bulk of the medium kpo- One can see 
that the potential changes much faster ai z < Zq than at z > Zq. The electric field can be 
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easily obtained now by differentiating Eq. (10), 



E{z) 



( Z — Zn . 



It is obvious that the maximum of the electric field is reached &X, z — Zq: 

Emax = Ci— - = Ci = Ci . (13) 

eA TT TT 

The areal charge density at 2; = 2^0 is given by eq. (9). Assuming that this charge is 
distributed over the layer of width A we can estimate the average charge density in this layer 



A TT^ 



epch ~ ^ = , (14) 



i. e. about 70% of the proton density at z — )■ —00. The energy density of the electric field 

is Sem — and therefore the total electrostatic energy is easily calculated as 

00 

E.^^ J S^^dV ^ ^ ci(2 - cO^Pp , (15) 

—00 

where S is the surface area and Pp is the bulk proton density. Note that this expression is 
approximately equal to the energy stored in a planar capacitor with areal charge density a+ 
and gap A. 

The authors of ref. [7] discuss the possibility of a new family of stable star-like objects 
(nuclear cores) where the repulsive electromagnetic force, associated with the charged layer 
at a sharp boundary, is balanced by the gravitational force. Here we present our counter- 
arguments concerning such possibility. First of all we point out that in a macroscopic nuclear 
matter, besides neutrons and protons, the electrons must be present to neutralize the bulk 
proton charge. By this reason the /5-equilibrium will be shifted to the neutron-rich side, 
so that the matter will contain mainly neutrons with a small admixture of protons and 
electrons. According to standard calculations of /3-equilibrated npe-matter in neutron stars 
(see e.g. refs. [11, 13]), the fraction of protons, ^ = Pp/ Pb should be on the level of a 
few percent. As well known, such matter is unbound, in contrast to the symmetric nuclear 
matter with Nn ~ Np without electrons. 

For simplicity we replace the baryon, proton and neutron densities by their mean values: 

PB = , Pp=^=^PB , Pr^^^^{l- QpB , (16) 



where R is the star radius. 

Now let us consider a star-hke configuration with sharp boundary at r — R, where we 
expect formation of a strong electric field. In our qualitative discussion we represent the 
star energy per baryon as the sum of three contributions: 

WiR) = W^UR) + WUR) + W^grav(i?) , (17) 

where terms in the r.h.s. represent the kinetic energy of degenerate fermions (npe), the 
electrostatic energy and the gravitational energy. The kinetic energy term is calculated 

easily for any fermionic species of mass rrii and Fermi momentum kpf 



(2x^ + Xi) ^1 + xf — arcsinh {xi] 



mi(l + ^xf), Xi^l ^^^^ 

1 



where i=n,p,e and Xi — kpi/mi. Below we include in VFkin only neutron contribution and 
express as a function of R using Eq. (16), 



The electrostatic energy of a spherical dipole layer can be obtained from Eq. (15) by 
replacing pp — ?,^Nb/^'kR^ and S — AnR^. Then the energy per baryon is 

WUR) = ^ci(2 - ^ , (20) 

where C3 ~ 0.426. It should be emphasized that this energy is by factor ^ ~ smaller 
than the energy of net charge Q — AirR'^a^ distributed in a spherical layer of wdth A and 
radius R (see ref. [8]). Finally, the gravitational energy of a baryon at the surface of a star 
of mass M — Nsms and radius it! is easily calculated in Newtonian gravity: 

GMniB f tubV Nb 



W = --^ = -^;^j (21) 

where mpi = 1/VG 1.22 • 10^^ GeV is the Planck mass and uib — 939 MeV is the nucleon 
mass. 

Putting together all the contributions we can now rewrite Eq. (17) as 



C3- - Nb 

e \mpij 
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where the ultrarelativistic hmit of Eq. (18) has been used in the kinetic term. One can 
immediately see that the electrostatic and gravitational contributions become equal at 



where the value ^ = 1/38 was used for the estimation. The corresponding mass Mmin = 
NrauJnB is about 10^^ g, i. e. almost 20 orders of magnitude smaller than the maximum 
mass of a neutron star. At Nb < -/Vmin there exists no equilibrium configurations because 
the kinetic pressure will force the system to expand indefinitely. 

At Nb > A^min the negative gravitational contribution exceeds the electrostatic one. 
The latter becomes completely negligible and can be neglected at baryon numbers ~ 10^^ 
characteristic for compact stars. As was first demonstrated by Chandraseckhar [22] and 
Landau [23], a compact star can be stabilized by the pressure of a degenerate fermion gas 
alone. Prom (22) one can see that the balance between kinetic and gravitational energy can 
be achieved if the total number of baryons is smaller than 



As follows from this simple calculations, the corresponding maximum mass of a neutron 
star, Mmax = A^max?7T.n IS about 3.2 Mq. Numerical calculations within General Relativity 
[24] for the degenerate noninter acting neutron gas predict somewhat lower maximum mass, 
Mjaax = 0.7Mq. The repulsive neutron- neutron interaction may increase this value up to 
4Mq (sec e.g. rcf.[25]). It is easy to sec that stars with A^min < Nb < A^max have a stable 
equilibrium state. Indeed, Eq. (22) shows that in this case W{R) > and changes as 1/R, 
so that the star should expand to lower densities. But this expansion will terminate at a 
large enough R, when the neutron gas becomes non-relativistic. As follows from Eq. (18), 
in this limit the energy per particle is m„ + 3A;|.^/10m„, i. e. changes as 1/R^ at large R. 
Since the negative gravitational contribution still changes as 1/R, it wins so that the total 
energy per baryon approaches m„ from below. Therefore a minimum develops in W{R). 
This is famous Landau's argument presented in ref. [23]. In ref. [25] it was reformulated for 
the case of General Relativity and for an arbitrary fermion mass. 

It is clear now that ignoring the kinetic-energy term in Eq. (22), as done in ref. [7], 
will lead to the erroneous conclusion regarding the star stabihty. Although the electrostatic 
and gravitational forces can balance each other at Nb — A^min, this state corresponds to 




(23) 




(24) 
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a maximum of energy. As one can easily see, the system will either expand or collapse 
depending on its actual baryon number compared to A^min- Moreover, the electrostatic force 
can support only the surface layer of matter where the electric field is nonzero. The inner 
layers of the star require the kinetic pressure to resist gravity. Therefore, the mechanism 
discussed in ref. [7] cannot produce any new stable object. Inclusion of the electrostatic 
effects will only slightly increase the star mass, see e. g. ref. [8]. 

In conclusion, we have demonstrated that very strong electric fields can be generated on 
a sharp boundary separating the bulk matter from the vacuum. The strength of the field is 
determined by the separation of positive and negative charges in the surface layer. For the 
proton charge distribution with the diffuseness parameter of a nuclear scale the field strength 
may exceed significantly the critical value for the electron-positron pair production. The 
field strength diminishes when the charge distribution becomes more and more smooth. We 
have also considered star-like configurations with a strong electric field on the boundary. 
Our simple analysis shows that such electric field alone cannot balance the gravitational 
force to produce stable objects. The stable compact star configurations appear only when 
the counter-pressure of degenerate fermions is taken into consideration. For a more accurate 
description of compact stars with strong electric fields one should perform calculations within 
General Relativity and include effects of strong interactions in the energy of the nuclear core. 
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